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Abstract 

Samples in the system Lu2-xY_^Si207 (0^x^2) have been synthesized following the sol-gel method and calcined to 1300 °C, a 
temperature at which the j^-polymorph is known to be the stable phase for the end-members Lu2Si207 and ¥281207. The XRD 
patterns of all the compositions studied are compatible with the structure of the j^-polymorph. Unit cell parameters are calculated as 
a function of composition from XRD patterns. They show a linear change with increasing Y content, which indicates a solid 
solubility of P-Y 2^1201 in f^-Lu2^\20i at 1300 °C. ^"^Si MAS NMR spectra of the different members of the system agree with the 
XRD results, showing a Hnear decrease of the ^^Si chemical shift with increasing Y content. Finally, a correlation reported in the 
Hterature to predict ^^Si chemical shifts in silicates is appHed here to obtain the theoretical variation in ^^Si chemical shift values in 
the system Lu2Si207-Y2Si207 and the results compare favorably with the values obtained experimentally. 
© 2004 Elsevier Inc. All rights reserved. 

Keywords: Solid solution; Y2Si207; Lu2Si207; ^^Si MAS-NMR; XRD; Group electronegativity; Polymorphism; Silicates 



1. Introduction 

Rare earth {RE) pyro silicates {RE2Si2^i) are well 
known for their polymorphism. They exhibit up to eight 
structure types (at normal pressure) depending on the 
temperature and the ionic radius of the RE cation [1,2]. 
All of them contain (Si207)^~ structural units (two 
silicon tetrahedra sharing one corner), except the a- 
polymorph which has the average pyrosilicate composi- 
tion but contains (Si04)^~ and (Si30io)^~ anions in 1:1 
ratio, i.e., isolated tetrahedra and chains of three Si 
tetrahedra. 

Many structural studies have been conducted on RE 
pyrosiHcates containing one type of RE cation. The 
investigation of the solid solubihty between the RE2 
Si207 structures containing different RE cations has 
been, however, very scarce [2]. The study of this subject 
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is interesting, not only from the academic point of view, 
but also in the materials science field. Mixtures of 
refractory RE oxides are used as sintering aids for the 
fabrication of silicon nitride structural ceramics. Upon 
anneahng, {RE)2^i20i crystallizes at the grain bound- 
aries of Si3N4, and improves the high temperature 
mechanical properties of the material [3-8]. Knowledge 
of the crystalhne structures adopted by the {RE)2^i20i 
intergranular phase at different temperatures and RE 
contents is therefore of great value in understanding the 
behavior of these materials. Several researchers have 
shown that the high temperature strength and oxidation 
resistance of Si3N4 are correlated with the radius of the 
RE cation in the oxide additives, such that the smaller 
the RE the better the properties of the silicon nitride 
[8,9]. Choi and Lee [9] showed that Y2O3, LU2O3, and 
SC2O3, among all lanthanide oxides, impart the 
highest flexural strength values to Si3N4. In this paper, 
we have selected the first two oxides to examine the 
phase transitions with composition in the system 
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Lu2Si207-Y2Si207; members of this system are expected 
to form in the intergranular region of Si3N4 when 
mixtures of both oxides are used as sintering aids of the 
Si3N4 structural ceramic. 

Y2Si207 shows up to five polymorphs with increasing 
temperature (y, a, ^, y and also called y, B, C, D and 
E, respectively) [2]; Ito and Johnson [10] established 
the following transition temperatures in Y2Si207: 
a — > P — > y <^ d. They did not include the other 
low temperature polymorph (y). Lu2Si207 exhibits, 
however, a unique polymorph (P or C) up to the melting 
point of the compound. The study of the Lu2_xY;^Si207 
system is therefore of special interest to analyze the 
stability fields of the different polymorphs with tem- 
perature and composition as well as to study the 
structural variations across the regions of soHd solubi- 
lity. In this paper, we have selected a temperature 
(1300 °C) at which both Lu2Si207 and Y2Si207 show the 
P polymorph as the thermodynamically stable phase, 
and we have studied the long- and short-range order 
variations of the structures appeared with increasing Y 
content from pure j5-Lu2Si207 to pure ^-Y2Si207. 



2. Experimental section 

2.1. Synthesis 

The sol-gel route used for this study was derived from 
the synthesis of a well-homogenized gel of Y2Si207 [11]. 
The starting materials were Y(N03)3 • 6H2O (99.9% 
Sigma), Lu(N03)3 • 6H2O (99.9% Sigma), Si(OC2H5)4 
(TEOS, 98% solution Sigma), HCl 35% aqueous 
solution and 96% ethanol. A TEOS solution in ethanol 
(1:3 in volume) was added over appropriate amounts of 
Y(N03)3 . 6H2O, Lu(N03)3 . 6H2O and HCl for the 
preparation of Lu2_xY;^Si207 members with x = 0, 0.2, 
0.5, 1.0, 1.5, 1.8 and 2.0. The mixtures were stirred at 
40 °C for 72 h and the transparent gels obtained were 
dried at 60 °C for 24 h in air. Nitrate was eliminated by 
calcination at 500 °C for 1 h at a heating rate of 



l°Cmin~^ The powder obtained in each case was 
subsequently calcined at a heating rate of 5°Cmin~^ 
upto 1300°C for 24 h and slowly cooled down to room 
temperature. 

2.2. Characterization 

The global composition of the samples was examined 
by X-ray fluorescence (X Siemens SRS-3000); Table 1 
shows the agreement between the nominal and real 
compositions. Several complementary techniques were 
selected which allowed a detailed ex situ characteriza- 
tion of the samples. X-ray diffraction (XRD) studies 
were carried out using a Siemens D-501 diffractometer, 
with Ni-filtered Cu^a radiation, steps of 0.02° and 
counting time of 10 s. The XRD patterns were analyzed 
using both the Le Bail and the Rietveld methods with 
the GSAS software (Larson and Von Dreele) [12]. The 
background was previously subtracted automatically 
using the diffractometer software although the back- 
ground intensities were subsequently refined. Lattice 
and profile parameters (linewidths and asymmetry 
parameters) were calculated by means of a Le Bail 
refinement. A structural refinement was then applied 
following the Rietveld method; the refined parameters 
were: background coefficients, phase fractions, lattice 
constants, site occupation fractions for the RE site and 
isotropic temperature factors. Powder elemental silicon 
was used as internal standard. ^^Si Magic Angle 
Spinning Nuclear Magnetic Resonance (MAS-NMR) 
spectroscopy measurements were carried out in a Bruker 
DRX400 (9.39 T) spectrometer equipped with a multi- 
nuclear probe, using 4 mm zirconia rotors spinning at 
1 1 kHz. A single pulse sequence was used, with a pulse 
length of 2.5 |is {n/2 pulse length = 7.5 |is), an observa- 
tion frequency for ^^Si of 79.49 MHz, and an optimized 
delay time of 600 s. The chemical shifts are reported in 
ppm from tetramethylsilane (TMS). The spectra were 
simulated using a modified version of the Bruker Winfit 
program to handle the finite spinning speed in MAS 
experiments [13]. 



Table 1 

FRX analyses, ^^Si chemical shifts, full-width at half-maximum and fraction of Gaussian line in the Gaussian-lorentzian function (from the 
simulation of the experimental MAS-NMR spectra in Fig. 4), group electronegativity sums for the members of the solid solution, predicted ^^Si 
NMR chemical shifts (following Janes and Oldfield correlation [19]) and difference between the calculated and observed values 

FRX d (obs) (ppm) FWHM (Hz) x [xG/(l - x)L\ EEN d (cal) (ppm) Error (ppm) 

X [Lu2_jcYxSi207] (nominal) x \Lvi2-^ x^^i^i] (real) 



0.00 0.00 -91.70 

0.20 0.23(1) -91.87 

0.50 0.56(1) -92.05 

1.00 0.97(1) -92.62 

1.50 1.51(1) -93.03 

1.80 1.79(1) -93.38 

2.00 2.00 -93.65 



117.6 0.05 15.25 -91.85 -0.15 

130.4 0.09 15.26 -92.00 -0.13 

152.6 0.32 15.27 -92.22 -0.17 

158.6 0.42 15.28 -92.58 0.04 

147.1 0.43 15.29 -92.95 0.08 

118.6 0.25 15.30 -93.16 0.22 

47.31 0.34 15.31 -93.31 0.34 
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3. Results and discussion 

Fig. 1 exhibits representative portions (15-50° 20) of 
the XRD diagrams of different \.U2- x^^i^i members 
(x = 0.0, 0.2, 0.5, 1.0, 1.5, 1.8 and 2.0) calcined at 
1300 °C for 24 h. The diagram of the x = 0.0 sample 
matches the standard pattern of ^-Lu2Si207 (JCPDS 
chart 35-326), as expected. Increasing Y content 
produces very similar patterns, with variations in peak 
positions and intensity; see, for example, the shift of the 
peaks under the dashed lines, the intensity decrease of 
the diffraction at ca. 20° 20 and the progressive 
separation of the double reflection at ca. 42.5° 20, 
marked with arrows in the figure. These variations result 
from changes in unit cell dimensions and composition. 
The diagram of the x = 2.0 \.U2- x^^i^i sample 
matches the standard pattern of j^-Y2Si207 (JCPDS 
chart 38-0440). The structure of the ^-polymorph is, 
therefore, maintained across the Lu2Si207-Y2Si207 
system at 1300 °C. Traces of ^2-RE2SiOs appear as a 
secondary product of the synthesis (corresponding 
reflections marked with asterisks). Fig. 2a shows the 
FWHM of several reflections for all the compositions. 
The crystallinity of the samples varies with the amount 
of yttrium, increasing from x = 0.0 to 0.5 (FWHM 
drops about 14%) and staying roughly constant up to 
the yttrium-rich end. 

In order to calculate the dimensions of the unit cell of 
each member of the system, the patterns, excluding the 
reflections of X2-i^£'2Si05, have been analyzed using 
the Le Bail and the Rietveld methods, as described in the 
experimental section. The starting parameters for the 
refinement of the structures with x < 1 have been taken 
from those reported for pure ^-Lu2Si207 [14] while the 
refinement of compositions with x > 1 was carried out 
from the parameters reported for pure ^-Y2Si207 [15]. 
Fig. 2b exhibits the observed pattern, the GSAS fitted 
curve and the difference curve obtained from the 
structural refinement of the x = 0.5 sample (the rest of 
compositions show very similar fits to this one). It can 
be observed that all the reflections could be fitted on the 
basis of a monoclinic unit cell with space group C2/m, 
as corresponds to the crystal structure of the ^- 
polymorph. The lattice parameters as a function of Y 
content are shown in Fig. 3. The same size of the j-axis 
has been used in the three plots to appreciate relative 
changes in the a, b and c unit cell parameters. Both a 
and b lattice parameters increase linearly, at approxi- 
mately the same rate, with increasing Y content, while 
the c parameter does not change appreciably with 
composition. A possible explanation for the constant 
value of the c unit cell parameter versus the increase of a 
and b can be understood from the crystal structure of ^- 
Lu2Si207 [14] represented in Fig. 3. It can be observed 
that the RE atoms (Lu in this case) are aUgned along the 
c-axis, with a Lu-Lu distance of 4.71 A. The O 




Fig. 1. Representative portions of the XRD diagrams of samples in 
the Lu2-xY;,Si207 system with x = 0.00, 0.20, 0.50, 1.00, 1.50, 1.80 and 
2.00, calcined at 1300 °C for 24 h. Three lines have been drawn on top 
of three diffraction Hnes of the Lu2Si207 XRD diagram in order to 
appreciate their shift to lower angles with increasing Y content. The 
arrow shows a double peak with separation increasing with increasing 
Y content. * = traces of^ X2-RE2SiOs. 



neighbors of the Lu atoms are, however, distributed 
along both, the a and b directions, with much shorter 
distances (2.20 and 2.21 A — irregular coordination 
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Fig. 2. (a) FWHM of reflections indexed in Fig. 1. (b) Experimental 
(crosses) and calculated (solid line) XRD patterns of the x = 0.5 
Lu2-xY;cSi207 sample. The difference plot is also shown. Circles show 
the diffractions from the siHcon standard. 



octahedra). It seems reasonable that, on replacing Lu 
by a bigger atom like Y, the c unit cell dimension 
does not experience any size modification because there 
is enough room along the c-axis for the Y atoms to 
accommodate. However, the size of the a and b 
parameters increase because the bigger Y atom expands 
the Y-O octahedron in the a and b directions of the 
unit cell. 

Finally, the ^ angle of the unit cell (Fig. 3) decreases 
Hnearly with increasing Y content while the unit cell 
volume increases. A linear increase in the unit cell 
volume with increasing the RE^^ ionic radius has also 
been reported for pure ^-RE2^i20i members {RE=\.u, 
Yb, Tm, Er and Ho) [16]. In conclusion, the linear 
behavior of the unit cell parameters with composition 
is a clear indication of the soHd solubiHty between 
jg-Lu2Si207 and P-Y2^i20i at 1300 °C. 

The local environment of the Si nuclei in the members 
of the j5-Lu2Si207-j^-Y2Si207 system has been studied by 
means of ^^Si MAS-NMR spectroscopy (the crosses in 
Fig. 4 show the experimental ^^Si MAS-NMR spectra). 
All the spectra show a main ^^Si resonance in the ppm 



range characteristic of the environments of the 
[Si207]^~ units in the pyrosiHcate structure. Some 
spectra exhibit a small contribution at ca. —78 ppm 
(not shown), from the secondary product X2-i^£'2Si05. 
It can be easily observed from Fig. 4 that the maxima of 
the curves shift to lower frequency with increasing Y 
content. In order to obtain the chemical shift and the 
FWHM of each resonance, the spectra were simulated 
successfully (solid curve in Fig. 4) on the basis of a 
unique Gaussian-Lorentzian Hne ( = xG/(l-x)L, where 
G is the Gaussian function and L the Lorentzian 
function). Resulting parameters are given in Table 1. 

The ^^Si chemical shifts do vary with composition in a 
linear manner, as seen in Fig. 5a. The ^^Si chemical 
shifts of the end members are in agreement with the 
values reported in the Hterature for j5-Lu2Si207 [17] and 
j^-Y2Si207 [18]. The linear behavior of the chemical 
shifts with composition also support the XRD data in 
favor of the soHd solubiHty between ^-Lu2Si207 and 
i5-Y2Si207 at 1300 °C. 

Janes and Oldfield [19] proposed a ^^Si chemical shift 
correlation based on group electronegativity. The 
following section of this work is devoted to check that 
correlation in our 1^-Lu2^i2^i-1^-^ i^^i^i system. Janes 
and Oldfield first defined three types of silicon depend- 
ing on the presence of only (7-(type S silicon), a- and n- 
(type P silicon) bonded Hgands or siHcon coordinated to 
both S and P Hgand types (type M silicon). For type P 
siHcon sites, as is the case for our pyrosiHcate 
Lu2_xY;cSi207 structures, they give the following 
correlation: 



^si = -24.336 ZEN + 279.27, 



(1) 



where ZEN is the group electronegativity sum. In order 
to use the method, the first step is to estabHsh the group 
electronegativities for the various oxy-metal fragments 
to which Si is bonded in the pyrosiHcate structure. 

The group electronegativity sum of a given silicon 
site is determined, following Janes and Oldfield [19], as 
foHows: 



i:EN(e^) = lENf + (4 - ^) i:(ENnfZnf /Iznf), 



(2) 



where n is the number of total framework tetrahedra 
coordinated to a given silicon tetrahedron; in our case, 
n= \ because we are dealing with pyrosiHcate struc- 
tures. ENf and ENnf represent the group electronega- 
tivity values for framework and nonframework Hgands, 
respectively, and Znf is the formal charge of a nonframe- 
work cation. In order to account for variations in the 
bridging-oxygen bond angles, which are beHeved to 
influence the oxygen orbital electronegativities, Janes 
and Oldfield [19] derived the OSi group electronegativity 
as a function of bridging bond angle (SiOSi) using three 
siHca polymorphs (low quartz, low cristobaHte and 
coesite). Albeit for only three minerals, a Hnear 
correlation of 0.99997 is observed between the group 
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Fig. 3. Unit cell parameters {a, b, c, P and volume) plotted as a function of Y content in the system Lu2Si207-Y2Si207. The error bars are 
approximately the size of the symbols. The stars are the values found in the literature for the end-members j5-Lu2Si207 [14] and ^^-¥281207 [15]. 
Down-right: Crystal structure of j6-Lu2Si207 based on the description by Soetebier and Urland [14]. Polyhedra correspond to Si04 tetrahedra while 
the spheres represent Lu atoms. 



electronegativity and bridging bond angle as follows: 
EN(OSi) = (ZSiOSi/136.79) + 2.9235. (3) 

We have used this correlation to obtain the OSi group 
electronegativity in the members of our solid solution. 
The SiOSi bond angle is 180° for all the structures across 
the i5-Lu2Si207-j5-Y2Si207 join, and therefore, EN(OSi) 
has a constant value of 4.24 and that will be the value of 
lENf in Eq. (2). 

Janes and Oldfield [19] obtained group electronega- 
tivities of several oxy-metal Hgands but they did not 
show values for OLu or OY. We have calculated them 
empirically for the series of model compounds Hsted in 
Table 2, which include three compounds with OY 
groups and two with OLu groups attached to Si. We 
have done the calculations using Eqs. (1) and (3) (the 
latter when necessary) as well as the ^^Si chemical shifts 



reported in the Hterature for those compounds. The 
average values obtained in this way for OY (3.69) and 
OLu (3.67) have then been used for the calculation of 
the electronegativity sums of the Si site in our 
Lu2_xY;cSi207 system: 

I:EN(Q1) = EN(OSi) + 3i:(ENnfZnf/Iznf) = 4.24 

+ 3{[;cEN(OY) 3/6] + [(2 - x)EN(OLu)3/6]} 
= 15.25 + 0.03X. (4) 

The resultant group electronegativity sums are given 
in Table 1 , and they provide the basis for the prediction 
of the chemical shifts for the different members of the 
solid solution, which is obtained by substituting expres- 
sion (4) in Eq. (1): 

3s, = -24.336 ZEN + 279.27 = -91.85 - 0.73x. (5) 
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Fig. 4. Experimental (crosses) and simulated (solid line) ^^Si MAS- 
NMR spectra of samples in the Lu2_;cY;cSi207 system. 



The ^^Si chemical shifts obtained from Eq. (5) are 
shown in Table 1 and plotted versus x in Fig. 5a. It can 
be observed that the predicted values compare very 
favorably with our experimental values for ^^Si. The 
maximum difference between the calculated and the 
experimental chemical shift values is 0.34 ppm, which is 
of the order of the lowest differences presented in the 
work of Janes and Oldfield [19]. Although Janes and 
Oldfield present data for 99 Si sites in 51 different 
compounds, only 5 sites are of type, the majority (67) 
being sites. The present work reconfirms the good- 
ness of their correlation, which provides chemical shift 
values for pyrosilicates with a minimum error. 

Finally, the FWHM values obtained fitting the 
spectra to a unique Gaussian-Lorentzian curve and 
given in Table 1, have been plotted versus composition 
in Fig. 5b. If the degree of crystallinity of the samples 
were the only factors affecting the FWHM of the ^^Si 
MAS NMR peaks, then they would exhibit a similar 
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Fig. 5. (a) Squares: Experimental Si chemical shift values plotted as 
a function of T content in the system Lu2Si207-Y2Si207. Dashed Hne: 
^^Si chemical shift values predicted from the correlation of Janes and 
Oldfield [19] using electronegativity group values, (b) FWHM values of 
the ^^Si MAS-NMR spectra as a function of Y content in the system 

LU2Si207-Y2Si207. 



behavior to that of the diffraction peaks (Fig. 2a). This 
is not the case, although it could explain the differences 
in FWHM values between both end-members. We think 
that the distribution of Y and Lu in the crystal structure 
of the P-RE2Si2^i pyrosilicate plays an important role 
in the FWHM values of the NMR peaks. Each silicon 
nucleus in the P-RE2Si20i structure is surrounded by six 
RE atoms. A homogeneous distribution of Y and Lu 
atoms in the crystal structure would cause a unique type 
of Si environment, Si(6Lu) and Si (6Y), in both 
Lu2Si207 and Y2Si207 end-members, respectively, as 
well as in the x = \ Lu2-xY;cSi207 composition -Si(3Lu, 
3Y)-. On the contrary, a random distribution of Y and 
Lu in the crystal structure would give rise to seven 
different types of Si environments in any composition, 
namely. Si [(6— x)Lu, xY], with x = 0 to 6, the 
proportions of each depending on the Lu/Y ratio. This 
variety of environments would give rise to a maximum 
FWHM value in the spectrum of the x = \ 
Lu2-xY;cSi207 composition while the FWHM value 
would decrease when crossing to both end members. 
The latter seems, therefore, to be the case of our system; 
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Table 2 

^^Si chemical shifts, mean SiOSi angles and group electronegativities for a series of rare earth sihcates 



Group 


Compound 


^^Si chemical shift 


Mean SiOSi angle 


Group 


Average group 






(ppm) 


(deg) 


electronegativity 


electronegativity 


OY 


;6-Y2Si207 


-93.7 


180 


3.69 


3.69 


OY 


7-Y2Si207 


-92.7 


172 


3.70 




OY 


X2-Y2Si05 


-79.4 




3.68 




OLu 


;6-LU2Si207 


-91.7 


180 


3.67 


3.67 


OLu 


X2-Lu2Si05 


-77.1 




3.66 





we have tried to fit each spectrum to the sum of those 
seven different Si environments but the quahty of the 
spectra, together with the small range of ^^Si chemical 
shifts expected for the seven environments (from —91.7 
to —93.7 ppm) did not allow to reach any rehable 
results. 



4. Conclusions 

A soHd solution between Lu2Si207 and Y2Si207 is 
formed at 1300 °C. Yttrium is able to replace Lutetium 
in the RE crystallographic site while maintaining the 
structure of the j5-polymorph. 

The a and b unit cell parameters increase linearly with 
increasing Y content from Lu2Si207 to Y2Si207, while 
the c unit cell parameter does not change with 
composition. Increasing Y content produces, as well, a 
linear increase in the unit cell volume and a linear 
decrease in the P monoclinic angle. The linearity 
observed for the crystallographic parameters indicates 
soHd solubility in the system j^-Lu2Si207-^-Y2Si207. 

The ^^Si chemical shift values increase linearly with 
increasing Y content. The experimental values coincide, 
within a minimum deviation, with those predicted from 
electronegativity group values with the correlation 
proposed by Janes and Oldfield [19]. The analysis of 
the FWHM of the spectra indicates that Y and Lu 
atoms are randomly, rather than homogenously, dis- 
tributed in the crystal structure of the pyrosilicate. 
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